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Abstract 

£ — ■ \ We demonstrate a novel approach that allows the determination of very general 

classes of exactly solvable Hamiltonians via Bethe ansatz methods. This approach 
combines aspects of both the co-ordinate Bethe ansatz and algebraic Bethe ansatz. 
The eigenfunctions are formulated as factorisable operators acting on a suitable 
reference state. Yet, we require no prior knowledge of transfer matrices or con- 
served operators. By taking a variational form for the Hamiltonian and eigenstates 
we obtain general exact solvability conditions. The procedure is conducted in the 
framework of Hamiltonians describing the crossover between the low-temperature 
phenomena of superconductivity, in the Bardeen-Cooper-Schrieffer (BCS) theory, 
and Bose-Einstein condensation (BEC). 

cn ! 1 Introduction 

o _ 

Recent work suggests an adequate understanding of the nature of the BCS-BEC crossover 
requires the study of Hamiltonians exhibiting coupling between Cooper pairs of atoms 
and bosonic molecular modes [H[2]. Further, the necessity of exact solutions to properly 
explain the physics of superconducting devices such as ultrasmall superconducting grains 
[3] motivates the search for exactly solvable generalisations of the BCS theory. 

A standard technique in this field is the implementation of the Quantum Inverse Scat- 
tering Method (QISM) [3H6] to construct quantum Hamiltonians with multiple conserved 
operators, and then find their exact solution by Bethe ansatz methods. Typically, a solu- 
tion of the Yang-Baxter equation (YBE) [7H9] is used to construct a transfer matrix which 
generates the conserved operators of the Hamiltonian. The exact solution may then be 
obtained via Bethe ansatz techniques. We note that in principle the implementation of 
the co-ordinate Bethe ansatz [10] does not depend on any prior knowledge of an associated 
solution of the YBE, nor the conserved operators of the Hamiltonian that it generates. A 
more modern approach in the framework of the QISM is the algebraic Bethe ansatz [1H6]. 

The strategy demonstrated here enables the direct determination of general classes 
of exactly solvable Hamiltonians, with both a bosonic and Cooper pairing degrees of 
freedom, such that they model BCS-BEC crossover behaviour. By contrast with standard 
techniques, this gives a unified construction for classes of exactly solvable Hamiltonians 
with multiple free coupling parameters. By taking appropriate limits we recover seven 
exactly solvable subcases from the literature [TTHTT]. 
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In Section [2] a general Hamiltonian that models the BCS-BEC crossover behaviour 
is introduced. In Section [3] the method of obtaining exact solvability constraints for 
the variational Hamiltonian is outlined. Calculations that have been omitted for brevity 
appear in the literature [IS] . 

2 Variational BCS-BEC crossover Hamiltonian 

We consider the general Hermitian family of pairing Hamiltonians coupled to a bosonic 
degree of freedom 

H = H -H 1 (1) 

where 

L 

H = aN + kNI + f(zk)N k , (2) 
k=i 

L L L 

fc=l fc=l k,s 

for some complex- valued function f(z) and real- valued function g(z), and real- valued 
parameters a, k, (3 and a, which will be subject to certain solvability constraints yet to 
be determined. An overline is used to denote complex conjugation. The operators b\ and 
N k = b[b k for k > are hard-core Cooper pair creation and number operators^. There 
is a single bosonic mode with operators b^ and N Q . The particle operators satisfy the 
following commutation relations: 

MJ = p},iJ] = o v,-,*>o, M'] = { {I J Nt)iik J 'l:° 

where 5j k denotes the standard Kronecker delta. The Hamiltonian commutes with the 
total number operator TV = No + ^2 k N k . Hamiltonians in this family describe a system of 
bosonic molecules, condensed into a single bosonic degree of freedom, coupled to L Cooper 
pairs which are bosonic-like pairs of fermions that must observe an exclusion principle. 
The Hamiltonian consists of a diagonal part H , given in Eq. [21 describing the allowed 
Cooper pair energy levels and the self-interacting bosonic mode, and a cross-interaction 
part Hi, given in equation Eq. [31 describing level dependent molecule-pair coupling and 
pair-pair couplings. For the special case a = (3 = k = 0, which suppresses any action 
of the Hamiltonian on the bosonic part of the underlying Hilbert space, the Hamiltonian 
Eq. [I] reduces to a general form of the BCS Hamiltonian. It is for this reason we can refer 
to the general pairing Hamiltonian as a BCS-BEC crossover Hamiltonian. 

We will attempt to directly solve the Hamiltonian by formulating the wave- functions as 
factorisable operators acting on a suitable reference state, analogous to the algebraic Bethe 
ansatz [1]. On the other hand, we note that a co-ordinate Bethe ansatz approach [T0lll9il20j 
need not resort to any prior knowledge of a transfer matrix or a set of conserved opera- 
tors. We do not expect that the Hamiltonian has an exact solution in general, however, 
by combining these aspects of the co-ordinate and algebraic Bethe ansatz methods, this 

1 The Cooper pair operators are related to the fermionic operators c\ via the identity b\ = c*_ k c k where 
fermions are paired such that there is zero total momentum. 
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technique allows us to obtain exactly solvable models in a very general fashion. In partic- 
ular, we have found solvability conditions for two sub-classes of the Hamiltonian, namely 
the cases i) k = of no self-interaction term and ii) a = of no BCS pair-pair scattering 
term. 

3 Variational Approach for Exact Solvability 

Motivated by the approach of Richardson [11], we assume the ansatz, 

M L 

l*> = II C ^)l >' C(y)= 1 (y)bl + Y,Ky,z k )bl (4) 

j=l k=l 

for the eigenstates of Eq. [TJ where |0) denotes the vacuum state, h(y,z) is yet to be 
determined, and y G C. In general, acting the Hamiltonian H on the state will lead 
to terms that are linearly independent of \^)\ 

H\m) = E\m) + \<&) , (tf|$) = 0. 

The objective is to choose appropriate constraints that allow these terms to be isolated. 
The ansatz will be an eigenstate if the coefficients of the terms linearly independent 
of 1^) cancel. The Hamiltonian is furnished with an exact solution on the manifold 
in the coupling parameter space obtained through compatibility of the corresponding 
constraints. Technical details of the calculation can be found in the literature [18]. We 
introduce the notation 

M M 

By direct calculation the action of H on the state is of the form 

ML ML 
j=l k=0 k,s=0 

for distributions Oi depending on the various functions introduced above. 

In order to determine the exact solution we look to solve the eigenvalue problem 
H\^>) = E\^/) for some scalar E. The set of solvability constraints required to find such a 
solution define a manifold in the coupling parameter space along which the Hamiltonian 
has an exact solution. We briefly remark on the specific choices for the constraints. 
Imposing a constraint of the form 

2 (x, y, w, z) = g(w)k(x, y)h(y, z) + g(w)k(y, x)h(x, z), 

for some function k(x,y) to be determined later, allows the reduction of linear combina- 
tions of b\b\\^ij) and frj^ol^u) terms in to linear combinations of and frjl^i) terms 
by use of the definition C(yi)\^fji) = To keep equations concise we have introduced 
the parameter z , however, we take h(y,z ) = j(y), g(z ) = 1. For k ^ ^ tr, this 
constraint yields a trivial exact solution due to the incompatibility of the coefficients of 
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b\.bl\^ij) and ^fcol^i) t erms - Thus the k = and a = cases are treated separately. A 
constraint is chosen such that the eigenvalues are of a standardised form, 

h(y, z) = p~\a(a - y) + /3 2 )^^, E = £ % , 

/(*) -2/ ^ 

Compatibility of these constraints defines k(x,y). A final constraint is obtained by im- 
posing that the remaining coefficients cancel. These are analogous to the Bethe ansatz 
equations and, subject to the compatibility of all constraints, the Hamiltonian is exactly 
solvable if the yj are their roots. 

For the case a = compatibility of the constraints leads to the constraining relations 
for exact solvability: 

Pg(zk)j(yj) 



L f(\ M 



yj - (a + «) + > — — = 2 > , c 2 = k. 

ti /&) - vo i£ v, - vi 

for arbitrary constants /3(or a), ci, and c 2 . 

In the case n = the solvability conditions are 



Cl a/3 = c 2 (aa+/3 )/3, c x = c 2 f(z k ) - g(z k )g(z k ), 
h{yj,Zk) = /3 [<j{a-y j ) + ft )—- J -, 



a(a-y,j) + P 2 f(z k )- yj Vj - y t ' 

for arbitrary constants (3(or a), c±, and c 2 . In either case the j(yj) will be fixed by any 
normalisation of \^f>). 



4 Final Remarks 

We have determined manifolds in the coupling parameters of Eq. [1] for which an exact 
solution exists. As shown in Fig. [H in appropriate limits of the general exact models 
found above eight subcases are recovered, seven of which are known [TiT - HT] . Remarkably, 
this method has enabled a unified construction for exact solutions of models that were 
previously independent. There is opportunity to extend the scope of this approach to 
other models [2TH23] by adopting an ansatz of a different form and relaxing the assumed 
separability of the pair-pair interaction. 
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BCS-BEC H5] BCS-BEC [15] BCS 




BCS-BEC [13] Dicke [12] BCS [U] BCS [17J BCS [IT] 

Figure 1: Hierarchy of exactly solvable models that have been recovered. Specific models 
can be obtained from the references, excluding three previously unknown cases. 
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